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1.0  INTRODUCTION 

This  report  describes  a  study  which  had  as  its  objective  the  evaluation 
of  pursuit  and  proportional  navigation  laws  for  use  in  closed  loop  simulations 
to  develop/evaluate  new  techniques  for  emitter  homing  [1] .  The  use  of  a 

pursuit  guidance  law  was  a  potential  candidate  since  only  stationary  targets 
are  considered.  For  targets  which  have  a  significant  velocity  proportional 
navigation  has  long  been  considered  the  acceptable  approach  [2,  3],  There  have 
been  other  guidance  laws  which  have  been  considered  including  modern  optimum 
control  schemes  [2,  3].  However,  these  two  laws  with  their  variations,  adap¬ 
tations  and  extensions  cover  a  broad  class  of  useage.  In  addition,  their  sim¬ 
plicity  (as  compared  with  many  "optimal"  systems)  allow  a  good  understanding 
of  their  structure  and  operation.  Furthermore,  many  of  the  modern  control  laws 
can  be  considered  as  "augmented"  proportional  navigation.  Finally,  proper  use 
of  these  laws  result  in  systems  which  can  admit  fairly  wide  variations  in  para¬ 
meters  without  serious  degradation  of  performance  from  the  nominal  design  con¬ 
ditions. 

One  of  the  important  aspects  of  this  study  is  that  the  sensor  used  as  the 
control  element  is  likely  to  be  a  strapdown  type,  that  is  it  is  fixed  to  the 
body  of  the  vehicle.  It  is  important  to  note  the  ramifications  that  this  has 
on  the  mechanization  of  the  control  system  and,  hence,  its  sensitivity  to  the 
motions  of  the  vehicle,  which  can  be  quite  large.  This  is  in  contrast  to  many 
situations  in  which  the  homing  sensor  is  gimbal  mounted,  isolating  it  from  body 
motion. 

For  this  study,  a  framework  has  been  selected  which  is  relatively  simple 
yet  includes  the  essential  features  of  the  two  control  laws.  Provision  has 
been  made  to  introduce  the  major  physical  parameters  of  a  homing  flight  vehicle 
including  its  inertial  and  aerodynamic  properties  as  well  as  a  representation 
of  the  stability  augmentation  features  of  the  autopilot,  i.e.,  rate  gyro  and 
accelerometer  feedback.  For  the  present,  the  response  has  been  limited  to  the 
linear  regime  for  the  final  10  seconds  of  flight,  in  which  the  velocity  and 
aerodynamic  parameters  can  be  reasonably  considered  constant  (except  for  the 
closing  range  between  the  vehicle  and  target) . 
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The  two  types  of  error  models  used  are  described  in  Section  2.0.  The 
first  is  a  deterministic  model  in  which  a  systematic  external  disturbance  or 
an  internal  instrumentation  error  is  introduced  and  the  resulting  response  of 
the  system  is  computed.  The  second  is  a  stochastic  model  in  which  the  external 
disturbances  are  characterized  by  their  spectrum  levels.  The  stochastic  model 
utilizes  an  adjoint  technique  [5]  in  which  the  sensitivity  of  the  system  to 
various  types  of  errors  are  computed. 


In  Section  3.0  numerical  results  of  six  types  of  vehicles  are  presented. 
(Note  that  the  use  of  type  includes  the  pitch  and  yaw  characteristics  of  the 
vehicle,  if  they  are  different.)  The  tables  of  error  coefficients  allow 
various  combinations  of  error  levels  to  be  studied  and  an  overall  miss- 
distance  in  feet  to  be  computed. 

Finally,  in  Section  4.0  some  conclusions  based  on  the  numerical  results 
are  presented. 
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2.0  ERROR  MODELS 

2.1  Deterministic 

The  coordinate  system  with  the  vehicle  is  shown  in  Fig.  2.1.  For  our 
analysis  small  angles  are  assumed,  that  is 

cos  Y  ~  1 

cos  0  ~  1 

cos  a  ~  1 

X  =  — ^ - 

where  y  =  flight  path  angle 

0  =  body  angle 

a  =  angle  of  attack 
A  =  line-of-sight  (LOS)  angle 

Since  we  will  be  considering  only  the  last  10  seconds  of  flight  it  is  also 
reasonable  to  assume  tnat  the  velocity  is  constant. 

The  equations  of  motion  of  the  vehicle  are 

-mvy  =  L 

JQ  =  M 

where  L  =  aerodynamic  lift 

m  =  mass  (=  W/g)  (slugs) 

W  =  weight  (lbs) 

2 

g  =  acceleration  of  gravity  (=  32.2  ft/sec  ) 

J  =  moment  of  inertia 

Q  =  0 

M  =  aerodynamic  moment 
The  lift  and  moment  can  be  expressed  as 

-  -  Z  <x  +  Z-6 
m  a  o 
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7  ■  V  +  M«6  +  V 


where 

Zct  =  CLaq  S/m 

Zt  =  CT  xq  S/m 
0  Lo 

Ma  ’  <*,’  SWJ 

M5  "  V  SWJ 

”q  '  V  ^Wi/J 

a  =  angle  of  attack 
q  =  l/2p  V2 

6  =  control  surface  rotation 

3 

p  =  air  density  slugs/ft  (=  0.00238  at  sea  level) 

2 

S  =  reference  area  (ft  ) 

l  =•  reference  length  (ft) 

C  =  lift  curve  slope  (per  radian) 

La 

C  =  lift  from  control  surface  deflection  (per  deg) 

Lo 

C  =  change  of  mement  with  angle  of  attack 
ma 

C  .  =  change  of  mement  with  control  surface  deflector 
mo 

C  =  change  of  mement  due  to  body  rotation 
mQ 

For  the  small  angle  assumption 
Z  =  -V  Y  =  -V(0  -  a) 

Using  the  definition  of  angle-of-attack  as 

«  '  t 

where  u  is  the  lateral  motion  of  vehicle  in  body  coordinates. 
Z  =  u  -  VQ 
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inhere  the  control  signal  is  the  filtered  LOS  rate. 

The  block  diagram  for  the  pursuit  case  is  shown  in  Fig.  2.2  and  for  the 
proportional  case  in  Fig.  2.3.  As  may  be  seen  there  are  four  error  inputs: 

1.  Wind  disturbance,  e  .  The  wind  disturbance  can  take  the  form 

of  any  distribution:  For  example,  for  a  wide  band  gaussian  fluc¬ 
tuation  the  filter  (bandwidth,  oO  will  produce  a  reasonable  re¬ 
presentation  of  the  wind  gust  if W the  bandwidth  is  properly  selected. 
With  a  wide  bandwidth  a  fairly  pure  wind  shear  disturbance  can  be 
introduced  into  the  system.  For  the  models  shown  in  Figs.  2.2  and 
2.3,  a  wind  shear  of  1  ft/sec/100  ft  will  be  used  to  generate  the 
basic  sensitivity  of  the  system  to  winds.  In  terms  of  the  time 
history,  the  wind  shear  can  be  expressed  as  =  0.01  Vt. 

2.  Gvro  error.  ep.  For  the  pursuit  navigation  case  the  gyro  error 
enters  only  through  the  feedback  path  of  the  autopilot.  On  the 
other  hand,  for  proportional  navigation  case  it  enters  not  only 
into  the  damping  feedback,  but  also  into  the  command  loop.  The 
rate  gyro  signal  is  integrated  to  obtain  a  measure  of  the  body 
attitude,  0  ,  which  is  subtracted  from  the  LOS  measurement  from 
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Proportional  navigation  model 


l^f^,  I  i^^^^vw^'Jtv-?  v*>  *“*'>  ^^iv^ta 


the  sensor  giving: 


(0  -  A)  -  6. 


-A  +  e. 


where  e^  is  the  difference  between  the  actual  body  position  and  that 
estimated  from  the  rate  gyro  signal.  Thus,  in  the  block  diagram, 

X-j  represents  the  error  0g.  For  the  basic  sensitivity  analysis  the 
error  will  be  computed  on  the  basis  of  a  gyro  error  of  l°/sec. 


Sensor  error.  e_.  The  sensor  error  can  be  a  noise,  drift  or  bias. 
The  noise  portion  of  the  error  will  be  considered  below  where  the 
stochastic  model  is  described.  The  deterministic  error  is  assumed 
to  be  a  bias  of  1°  for  establishing  the  basic  error  sensitivity 
for  deterministic  errors. 


Accelerometer  error,  e. .  The  accelerometer  error  enters  the  feed- 

-  •  .  «  ^  .  4  t**..  .  t  _ , _ J  1  *  1.  _  T-1 _  .  .  .  I*  .  t _ _ _ 


back  path  of  the  acceleration  autopilot  loop.  For  some  of  the  cases 
there  will  be  no  accelerometer,  therefore,  k  =  0  and  there  will  be 
no  acceleration  error.  The  noise  component  wiM  be  considered  below 
and  the  steady  state  level  will  be  a  bias  of  1  ft/sec^  (about  30 
millig's) . 


The  sensor  transfer  function  is  assumed  to  be 


H„(S)  =  - - y 

s  (~  +  ir 


for  the  pursuit  case  and 


H.(S)  =  ? 

s  (—  + 1  r 

0) 

s 


for  the  proportional  case,  where  S  is  the  Laplacian  operator. 


The  equations  represented  by  the  block  diagram  can  be  cast  into  a  matrix 
equation  as 


X  =  FX  +  Gc 


where  the  definitions  of  the  components  of  this  matrix  equation  are  given  in 


Tables  2.1  and  2.2  for  pursuit  and  proportional  cases,  respectively. 


TABLE  2.1 


SUMMARY  OF  PURSUIT  NAVIGATION  EQUATIONS 


x  =  Fx  +  Ge 


X1  =  2 
x2  =  z 

x3  =  u 

x4  “  Q 

X5  "  A1 

X6  =  X2 
x^  =  0 

X8  =  UW 


el  eW 
£2  =  £S 
£3  =  £A 
£4  =  £G 


G23  =  k  KAZ6A 
G24  "  KQZ6A 
G33  =  k  KAZ6A 
G34  =  KQZ6A 
G43  =  k  KAM6 
G44  =  KQM6 
G52  =  WS 

G81  =  “w 


F12  =  1 


F23  =  V  A’  F24  KQZ6Aj  F25  =  G2KAZ<$A’  F28  "  "F23 
Z 

F33=A  F34  =  KQZ6A  +  V’  F35  =  G2KAZ6A’  F38  "  -F33 

M  Z  k  KM 

F43  =  r  +  -V^’  F44  "  MQ  +  KQV’  F45  =  G2KaV’  F4 
F55  =  '  “S’  F56  =  WS 


F61  Vtx  ’  F66  "‘V  F67  WS 

F74  =  1 
F88  =  '“W 

A  =  1  +  Zgk  Ka,  =  tQ  -  t  tQ  =  flight  time 
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The  numerical  solution  of  these  equations  can  be  expressed  as 

x  =  (I  +  AtF  )  X  ,  +  At  G  e 
n  n  n-1  n 

where  I  is  an  identify  matrix  and  At  is  the  solution  time  interval. 

2.2  Stochastic 

The  sensitivity  of  the  miss  distance  to  stochastic  inputs  is  most  conven¬ 
iently  done  using  the  adjoint  technique  [5].  Basically,  this  technique  recasts 
the  block  diagram  of  the  basic  system  into  an  adjoint  form  by 

1.  replacing  t  by  t  -  t  in  the  arguments  of  all  time  varying  coef¬ 
ficients  where  t  is  the  final  flight  time 

2.  reversing  all  signal  flow,  redefining  branch  points  as  sum  points 
and  sum  points  as  branch  points,  making  the  original  inputs  appear 
as  outputs  in  the  adjoint  system. 

Before  going  on  to  the  more  complex  problem  we  will  consider  a  simple 

system  to  illustrate  the  adjoint  technique.  Figure  2.4  shows  a  second  order 

system  with  its  adjoint.  The  system  is  excited  by  a  white  noise  input  e  having 

2 

a  spectral  density  of  $  units  /Hz.  The  rms  response  to  the  basic  system  is 
given  as : 


-hi  i«“)i2d“ 


where  H(<j)  =  — ^ - 2 - 

u)  -  oj  +  2£w  ju) 
o  o 

classic  result 


.  Carrying  out  the  integration  we  obtain  the 


For  the  adjoint  problem  the  solution  is  expressed  as 


.f  4  2 

=  S’#  w  x, 
J  o  1 


t.jUnaw  »-CV* I  .*«.»  f  l*VV*£1.. 


-Cw  t 


where  x  =  i  e  0  sinBt,  B  =  w  /l-?^  for  initial  condition  of  x. (o)  =  1, 
■*■  p  o  1 

x^(o)  =  0  (or  the  impulse  response). 

Carrying  out  the  integration  we  get 


©»> 


U) 


j  =  <v  sin  gt  _  - 

X1  (2(1-?2)  2>/I-?2 


sinBt  cosBt 


-2?w  t 
o 


U)  -2?w  t 

-^(e  °  -1) 

which  is  the  adjoint  time  history  of  the  variance  of  x^.  It  is  readily  seen 
that  as  2?Wot  becomes  large,  we  obtain  the  steady  state  result. 

°xl  =  4?  $ 


The  usefullness  of  the  adjoint  technique  becomes  apparent  for  more  complex 
systems  with  multiple  stochastic  inputs.  Not  only  are  all  the  sensitivities 
to  the  inputs  determined  in  one  run,  but  the  stability  of  the  system  may  also 
be  observed. 


Although  the  adjoint  technique  can  handle  nonlinearities  [5],  only  linear 
cases  will  be  considered  in  the  present  study. 


Using  the  two  rules  for  obtaining  an  adjoint  block  diagram  for  a  specific 
system,  the  pursuit  and  proportional  navigation  were  obtained  as  shown  in  Figs. 
2.5  and  2.6.  From  these  block  diagrams,  the  equations  for  both  systems  were 
obtained  in  the  form 

x  =  FAx 


where  the  x  is  a  9x1  vector.  The  r  is  the  system  adjoint  matrix  with  dimen¬ 
sions  appropriate  to  the  x  vector.  The  second  equation  is  a  vector  equation 
for  the  outputs  in  terms  of  the  solution  vector,  x.  It  is  of  4x1  dimension  to 
account  for  the  four  output  variances  (stochastic  inputs  for  the  original  sys- 
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TABLE  2. A 

SUMMARY  OF  PROPORTIONAL  NAVIGATION  EQUATIONS 


ADJOINT  SYSTEM 


x  =  F^x  S  =  SQ2 


=  0 

S„, 

=  (0  x, 

w  w 

01 

w  7 

■  “sH 

S02 

=  0)s2x5  -  CO, 

■  °a2/\ 

S03 

=  z6kKAA(x2 

-  “oH 

S04 

=  x9/ 

•  -  -V  „  A 

Vt  *16 

w 

s 

-  Vt 

tg 

=  1 

V  p  A 

v  A  ,  F33  - 

F  A 
*32  ’ 

f3a-^  + 

34  v 

F42-z5KqA>  F4'V  +  F42-  F4t-M6V  +  MQ 


42  6  Q  ’  43  “ 

17  ~  v  r>  A  -r»  A- 


F52  S  W>  F53  =  F52>  F54  =  M6KAG1A’  F55  =  ^>S 
F65  =  “s  ’  F66  =  *<V  F75  =  "^S  *  *76*  WS 


/6~  WS 


p  A  =  _p  A  p  A  =  p  A  -p  A  _  -F  A  p  A  _ 

82  *32’  *83  *82 ’  *84  *34’  *88  _ 


F99  ’  -l'V  F92  *  KqV'  F9t  *  F-  F93  *  F92’  F94  =  KqV 


A  =  1  +  kKAz„ 
A  o 


iiBWRWiigwsew^ 


tern).  The  details  of  the  equations  are  given  in  Tables  2.3  and  2.4  for 
pursuit  and  proportional,  respectively.  It  is  noted  that  a  filter  has  been 
added  between  the  gyro  noise  output  and  the  system.  This  is  to  account  for 
the  reasonable  long  correlation  times  that  can  exist  for  gyro  noise. 

Solution  of  these  equations  is  carried  out  by  seeking  the  impulsive 
condition  on  miss  distance  which  can  be  accomplished  by  giving  an  initial 
condition  of  z(o)  =  1  and  propagating  the  results  for  the  solution  time,  which 
was  taken  as  10  seconds  for  all  cases  run.  As  for  the  basic  system  described 
above,  these  equations  were  solved  numerically  as 


x  =  (I  -  AtF. )x  . 
n  A  n-1 


S  =  S  +  (S  2  +  S  2) 
n  n-1  2  o  o' 

n  n-1 


with  At  =  0.01  second  and  initial  condition  X2(o)  =1.  To  obtain  variance  of 
the  miss  distance  from  all  four  error  sources  we  multiply  each  output  by  its 
spectral  density  as 

a2  =  S.4>  +  S.4>  +  S,«.  +  S.$„ 
lw  2  S  3A  4G 

where  the  subscripts  are  wind,  sensor,  accelerometer,  and  gyro, 
respectively.  In  the  case  of  the  wind  spectrum  we  have 


$  =  — 
w  0) 


where  aw  is  the  rms  level  of  wind  gust 

w  =  V/L 
w  w 

L,,  is  the  wind  correlation  length 
w 

V  is  the  vehicle  velocity. 

Thus,  an  error  budget  for  the  system  is  quite  readily  obtained. 
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3.0  NUMERICAL  RESULTS 

Calculations  for  both  pursuit  and  proportional  navigation  laws  were  carried 
out  for  six  representative  flight  models. 

MISSILE  -  is  a  boost-glide  vehicle  with  a  speed  less  than  Mach  1 
in  the  terminal  phase  of  the  flight.  Since  it  is  symmetrical, 
both  pitch  and  yaw  are  essentially  the  same.  The  characteristics 
are  similar  to  the  Maverick  missile. 

PROJECTILE  -  is  based  on  a  155  mm  cannon  launched  projectile  which 
deploys  a  set  of  fins  and  wings  for  control.  It  also  has  a  speed 
less  than  Mach  1  in  the  terminal  phase  of  flight.  As  for  the  mis¬ 
sile,  it  also  is  symmetrical. 

MINI  DRONE  -  is  a  representative  version  of  a  powered  mini  drone. 

The  vehicle  considered  in  this  study  is  a  skid-to-turn  vehicle 
with  side  force  surfaces  for  enhanced  yaw  steering  capacity.  Its 
pitch  and  yaw  dynamic  characteristics  are  substantially  different 
from  one  another  so  both  pitch  and  yaw  models  are  considered. 

GLIDE  BOMB  -  bears  a  representative  likeness  to  the  GBU-15  plana? 
wing  glide  bomb.  It  also  is  a  skid-to-turn  vehicle  but  with  no 
enhancing  side  force  surfaces.  The  pitch  and  yaw  characteristics 
are  sufficiently  different  from  one  another  so  that  each  is  con¬ 
sidered  as  a  separate  case  in  the  error  studies. 

Table  3.1  gives  the  numerical  parameters  [6]  for  the  six  cases  studied. 

The  yaw  parameters  have  been  modified  in  sign  to  allow  the  equations  developed 
in  Section  2.C  to  be  used.  The  flight  characteristics  of  each  vehicle  required 
stability  augmentation  to  be  adequate  for  the  application.  Representative 
values  of  feedback  gains  were  selected  to  represent  this  stability  augmentation 
or  autopilot  function: 

K„  is  a  rate  gyro  which  senses  the  body  rate  and  supplies  a  signal 

to  the  control  surface  to  oppose  body  motion  and,  thus,  provides 
additional  damping. 

k  is  an  accelerometer  which  senses  the  body  acceleration  aid  provides 
a  signal  to  the  control  surface  which  is  used  to  null  the 
difference  between  commanded  acceleration  and  actual  body 
acceleration.  For  the  Projectile  and  Mini  Drone,  there  is 
no  accelerometer  so  that  the  gain  is  taken  to  be  zero. 

is  the  gain  which  converts  the  commanded  acceleration  into  a 
control  surface  position 
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TABLE  3.1 


PARAMETERS  FOR  ERROR  ANALYSIS 


MISSILE 

PROJECTILE 

MINI 

(PITCH) 

DRONE 

(YAW) 

GLIDE 

(PITCH) 

BOMB 

(YAW) 

S  (FT2) 

0.786 

0.196 

18.18 

18.18 

16.6 

16.6 

W  (LBS)  2 

398 

135 

135 

135 

3000 

3000 

J  (SLUG  FT  ) 

53.5 

6.13 

25.5 

32.4 

646.2 

716.6 

V  (FT/SEC) 

700 

800 

200 

200 

600 

600' 

l  (FT) 

1 

0.5 

1.65 

1.65 

1.54 

11.33 

l(ft) 

c.  (/RAD) 

CT?  (/DEG) 

(/RAD) 
c  (/DEG> 

900 

900 

900 

900 

900 

900 

-14.32 

-20.05 

-5.5 

-1.5 

-5.73 

-1.66 

-0.066 

-0.130 

-0.0066 

-0.0044 

-0.008 

-0.015 

-6.99 

-19.19 

-1.72 

-0.80 

-0 

-0 

-0.230 

-0.390 

-0.0223 

-0.140 

-0.038 

-0.006 

(/RAD) 

MQ 

-573 

-155 

-7.0 

-10.0 

-0 

-0 

k 

-1 

0 

0 

0 

-1 

-1 

Ka  (D/FT/S2) 

2.2/ g 

2/g 

2/g 

2/g 

2/g 

2/g 

V,  (D/R/S) 

K^;  (D/R/S) 

0.2*r 

0.2*r 

0.2*r 

0.2*r 

0.2*r 

0.2*r 

0.3*' 
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The  dynamic  characteristics  of  the  vehicle  autopilot  were  computed  and 
are  also  included  in  Table  3.1.  As  may  be  seen,  there  is  a  wide  variation  in 
nt  .ural  frequency  from  a  high  of  about  16  radians/sec  to  a  low  of  about  2.6 
radians/sec,  or  a  factor  of  6.  It  is  to  be  noted  that  there  are  differences 
between  the  dynamic  characteristics  of  the  vehicles  in  the  pursuit  mode  from 
those  in  the  yaw  mode.  The  reason  for  this  is  that  for  an  adequate  behavior, 
the  rate  damping  had  to  be  greater  for  pursuit  than  for  proportional,  except 
for  the  mini  drone.  The  aerodynamic  gain  varies  from  0.276  to  0.77  at  about  a 
factor  of  three.  The  overall  kinematic  gain  of  the  systems  varies  with  time  as 


G  ,  x  G. 
oi  i 


'kinematic 


t  -  t 
o 


where  i  =  1  is  the  proportional  navigation  case 
i  =  2  is  the  pursuit  navigation  case 
=  G^/V  the  navigation  gain 
Goi  -  aerodynamic  gain 
tQ  =  flight  time 


Although  the  gains  selected  are  not  necessarily  optimal,  they  do  provide  a 
reasonably  good  system  response  for  no  errors  or  disturbances. 


The  parameters  of  Table  3.1  were  used  in  the  numerical  evaluation  of  the 
equations  outlined  in  the  previous  section.  As  a  baseline,  the  six  venicles 
were  run  for  both  proportional  and  pursuit  with  no  errors  or  disturbances. 
Figures  3.1,  3.2,  and  3.3  are  typical  responses  of  the  cross  range  distance 
(miss  distance)  and  lateral  acceleration  for  the  Missile,  Mini  Drone-pitch, 
and  Glide  Bomb-yaw,  respectively.  The  initial  condition  in  all  cases  was  a 
10  ft/sec  cross  range  velocity,  z(o),  or  an  angular  miss  alignment  of  £(o)/V 
for  the  vehicle.  Except  for  the  pursuit  case  of  Glide  Eomb-yaw,  miss  distances 
were  small  (as  typified  by  the  Missile  and  Mini  Drone).  Because  of  the  poor 
Glide  Eomb-yaw  response  characteristics,  it  is  not  too  surprising  that  diffi¬ 
culties  existed  in  this  final  10  second  flight  period.  It  is  further  to  be 
noted  that  the  peak  acceleration  levels  of  the  pursuit  cases  were  constantly 
higher  than  for  the  proportional  cases. 
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Fig.  3.1.  Generic  missile  response  for  no  errors  -  initial 
velocity  =  10  ft/sec. 
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Table  3.2  gives  a  complete  summary  of  the  error  sensitivities  of  the  six 
vehicle  types  for  both  proportional  and  pursuit  navigation.  The  two  sections 
of  the  table  are  the  basic  error  sensitivities  for  deterministic  and  stochastic 
inputs  as  discussed  in  the  previous  section.  For  the  deterministic  cases,  the 
numbers  represent  the  miss  distance  in  feet  due  to  a  specific  level  of  input  as 
follows 


SENSOR 

WINDSHEAR 

ACCELEROMETER 

GYRO 


-  1°  (BIAS) 

-  1  ft/sec/100  ft. 

-  1  ft/sec2  fRIAS) 

-  l°/sec  (BIAS) 


For  the  stochastic  cases,  the  values  in  the  table  represent  the  variance  divided 
by  the  spectrum  level  as  noted  in  the  previous  section. 

In  looking  at  the  error  sensitivities,  certain  patterns  emerge  for  the  two 
control  laws.  If  we  leave  out  the  Glide  Bomb-yaw  case  as  being  atypical,  we 
notice  that  both  accelerometer  and  gyro  errors  are  relatively  low  contributers 
for  both  pursuit  and  proportional  for  deterministic  inputs.  On  the  other  hand, 
the  wind  shear  and  sensor  inputs  are  substantial  contributers  for  pursuit  and 
very  small  contributers  for  proportional.  When  we  look  at  the  stochastic  in¬ 
puts,  we  see  that  the  major  contributers  for  proportional  are  the  sensor  and 
the  gyro.  The  reason  for  the  large  gyro  contribution  is  due  to  the  fact  that 
the  gyro  is  used  to  separate  out  the  body  position  from  the  line-of-sight .  The 
wind  and  accelerometer  contributions  are  relatively  small.  For  pursuit,  the 
sensor  is  a  moderate  contributer  while  the  wind  gustiness  is  the  major  contri- 
buter. 

Figure  3.4  illustrates  some  typical  adjoint  responses  for  the  gyro  noise 
error  and  proportional  navigation.  The  data  have  been  normalized  with  respect 
to  the  maximum  values  given  in  Table  3.2  (for  example,  159  ft  /spectrum  level 
for  the  missile).  This  figure  provides  a  graphic  illustration  of  the  sluggish¬ 
ness  of  the  Glide  Bomb  (particularly  yaw)  compared  with  the  other  vehicles. 

It  is  further  noted  that  the  values  for  gyro  noise  given  in  Table  3.2  as¬ 
sume  wide  band  input.  These  values  are  markedly  reduced  as  the  correlation 


27 


bsSSBSSBBBHBI 


time  increases. 


The  relationship  shown  in  Figure  3.5  can  be  expressed  as 


where 


(£  \  (^fi 

Vo)  '  i  +  t* 

/oh 

I  —  I  is  the  value 

\V 

/  OG  \  ,  t„ 

I J  is  the  corre 


is  the  value  for  the  gyro  noise  in  Table  3.2. 


is  the  correlated  value 


Tg  is  the  correlation  time  of  the  gyro  noise 

This  expression  can  be  used  for  all  of  the  vehicles  as  a  reasonable  approxima¬ 
tion.  The  spectrum  level  for  the  gyro  noise  is  given  as 


<j>  —  2  o  T 

G  G  G 

To  obtain  a  measure  of  overall  performance,  it  is  necessary  to  select  some 

reasonable  values  of  errors/disturbances.  For  the  sensor,  we  have  assumed  a 

2 

bias  of  1  degree  and  noise  level  of  0.08  deg  /Hz.  For  a  bandwidth  of  10  Hz  this 
implies  an  rms  value  of  about  1.6°.  The  wind  shear  was  taken  as  C.5  ft/sec/100 
ft.  which  implies  a  variation  of  150  ft/sec  from  sea  level  to  30,000  ft.  altitude 
(for  the  pitch  plane).  The  gust  intensity  was  chosen  to  be  10  ft/sec  (rms),  with 
the  spectrum  level  being 

2  L  (10) 2  200L 

.  w  w 


where  Lw  is  the  correlation  length  assumed  to  be  900  ft.  and  V  is  the  vehicle 

velocity.  [7]  suggests  that  an  autopilot  grade  accelerometer  will  have  a  bias 

2 

of  10  millig's  (or  0.32  ft/sec  ).  For  our  analysis,  we  will  assume  a  pessimistic 
2 

value  of  1  ft/sec  bias.  [8]  suggests  that  autopilot  grade  accelerometers  are 

an  order  of  magnitude  lower  in  accuracy  than  inertial  grade,  and  that  a  noise  spec- 

2  2 

trum  level  of  1  millig  (or  0.024  (ft/sec  )  /  Hz)  will  be  reasonabley  conservative. 
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NORMALIZED  GYRO  NOISE  ERROR 
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Fig.  3.5.  Normalized  gyro  noise  error  as  a  function 
of  correlation  time. 
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[8]  suggests  that  the  drift  or  bifts  of  an  autopilot  grade  gyro  will  be  two  or¬ 
ders  of  magnitude  worse  than  an  inertial  grade  unit.  Considering  a  10-15  deg/h 
inertial  grade  gyro  implies  a  bias  level  of  about  0.35  deg/h  for  an  autopilot 
grade  gyro,  which  is  in  agreement  with  the  range  of  autopilot  grade  gyros  sug¬ 
gested  by  [ 7] . 

Finally,  a  chart  in  [7]  suggests  that  for  an  inertial  gyro  an  rms  of  0.1°/h 
with  a  correlation  time  of  15  min.  would  be  a  reasonable  noise  level.  Using  a 
factor  of  100  on  each  of  these,  we  obtain  a  spectrum  level  of 

G  =  j  (10)  =  1.5  x  10"4  (deg/sec)2/Hz 

With  a  correlation  time  of  10  seconds,  the  reduction  factor  for  the  gyro  noise 

-2 

sensitivity  in  Table  3.2  will  be  10  .  These  levels  make  the  gyro  noise  contri¬ 

bution  negligible. 

The  error  budget  using  these  values  is  summarized  in  the  two  sections  of 
Table  3.3  for  proportional  and  pursuit,  respectively.  The  values  for  the  indi¬ 
vidual  components  of  the  error  budget  are  the  10  miss  distances  in  feet  and  the 
RSS  error  for  each  of  the  vehicles  for  each  law  is  shown  at  the  right  hand  side 
of  the  table.  It  can  be  seen  that  the  RSS  errors  for  proportional  are  consistent¬ 
ly  lower  than  for  pursuit  by  a  substantial  measure.  For  proportional,  the  major 
contributor  is  the  sensor,  while  for  pursuit,  the  major  contributers  are  the  sen¬ 
sor  bias,  wind  shear,  and  gust. 

Although  Table  3.3  does  not  seem  to  indicate  so,  gyro  bias  can  be  a  sub¬ 
stantial  contributor  to  miss  distance.  The  reason  for  this  is  demonstrated  in 
Figure  3.6,  which  is  a  plot  of  cross  range  distance  and  lateral  acceleration  for 
the  missile  with  a  bias  of  l°/sec  in  the  gyro.  The  rather  large  excursions  and 
acceleration  introduced  by  this  error  stand  out  in  comparison  to  the  zero  error 
case.  The  accelerations  caused  by  gyro  bias  will  produce  a  substantial  increase 
in  drag  which  will,  in  turn,  reduce  the  speed  of  an  unpowered  vehicle.  The  re¬ 
duced  speed  will  reduce  both  range  and  maneuver  capability,  thereby  increasing 
the  potential  to  miss  the  target. 
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In  order  to  obtain  a  measure  of  the  sensitivity  to  the  gain  level,  a  sec¬ 
ond  case  was  run  in  which  the  overall  navigation  gain  was  reduced  by  25%.  The 
results  of  this  calculation  are  shown  in  Table  3.4.  As  for  Table  3.2,  the  two 
sections  give  the  basic  sensitivities  for  pursuit  and  proportional,  respectively. 
In  Table  3.5,  the  two  sections  give  the  error  budget  to  the  specified  input 
values  as  before.  Comparing  the  zero  error  case  for  the  two  gain  levels,  it 
may  be  seen  that  the  only  substantial  change  was  for  the  GBU-yaw  case  for  pro¬ 
portional,  in  which  the  miss  distance  was  somewhat  increased. 

Looking  at  the  individual  error  sensitivities,  we  see  that  for  the  propor¬ 
tional  case  that  error  sensitivities  for  sensor  and  gyro  noise  were  substan¬ 
tially  decreased,  while  the  sensitivity  for  gyro  bias  was  noticeably  increased. 
The  acceleration  noise  sensitivity  was  also  somewhat  increased,  but  it  was 
quite  low  in  any  event.  For  the  pursuit  law,  there  was  a  slight  degradation 
for  all  contributors.  The  total  effect  was  a  slight  decrease  in  the  RSS  error 
for  the  proportional  case  and  a  slight  increase  in  the  error  for  pursuit  case. 
Thus,  we  may  conclude  that  both  laws  can  perform  adequately  for  substantial 
changes  in  gain  from  the  nominal. 


(1)  HISS  DISTANCE  IN  FT  (la)  TYPICAL 


4.0  CONCLUSIONS 


The  following  observations  can  be  made  from  the  numerical  results  of  this 
study: 

1.  The  peak  accelerations  are  generally  higher  for  the  pursuit 
law.  This  means  that  the  potential  for  saturation  of  the 
maneuver  capability  of  vehicle  is  more  likely.  In  addition, 
a  greater  velocity  loss  can  be  expected  from  the  pusuit  law. 

2.  The  largest  error  contributors  for  pursuit  navigation  are  sen¬ 
sor  bias,  wind  shear,  and  wind  gust.  On  the  other  hand,  the 
proportional  law  is  relatively  insensitive  to  wind  effects. 

3.  The  largest  error  contributors  for  the  proportional  law  are 
sensor  noise  and  gyro  bias.  The  contribution  of  gyro  bias 
arises  from  the  need  to  subtract  body  rate  from  the  strap- 
down  sensor  measurement.  The  larger  sensitivity  of  the  pro¬ 
portional  law  to  sensor  noise  is  basically  due  to  the  fact 
that  noise  is  differentiated.  On  the  other  hand,  for  this 
same  reason,  the  pursuit  law  is  relatively  insensitive  to 
sensor  and  gyro  noise. 

4.  The  Mini  Drone  vehicle  with  proportional  guidance  produced 
the  best  performance  with  respect  to  miss  distance  in  this 
study.  In  this  case,  even  sensor  and  gyro  contributions 
were  small.  The  basic  reason  for  this  is  that  the  Mini  Drone 
vehicle  is  highly  maneuverable.  On  the  other  hand,  its  slow 
speed  may  make  it  undesirable  for  some  applications. 

5.  A  reduction  in  the  gains  of  both  proportional  and  pursuit 
laws  resulted  in  some  degradation  of  the  response.  For 
proportional  law,  the  gain  reduction  reduced  the  sensitivity 
to  gyro  and  sensor  noise,  but  increased  the  sensitivity  to 
gyro  bias. 

These  obser”'  ons  support  the  following  conclusions: 

-  !.>■*  accurate  results  in  closed  loop  simulations  of  systems 

with  proportional  navigation,  the  simulations  should  include 
the  effects  of  gyro  and  sensor  noise.  For  simulations  in¬ 
volving  pursuit  law,  accurate  simulations  should  include  the 
effects  of  the  wind  (both  shear  and  gust)  as  well  as  any  bias 
in  the  sensor. 


5®SS 


The  numerical  results  and  the  above  associated  observations 
support  the  conclusion  that  a  proportional  based  law  is 
generally  more  effective  than  a  pursuit  based  law. 


In  the  implementation  of  a  strapdown  sensor  for  proportional 
navigation  a  gyro  of  less  noise  and  drift  than  the  usual  auto¬ 
pilot  sensors  may  be  needed. 


Closed  loop  simulations  [6]  of  these  vehicles  with  a  more  accurate  repre¬ 
sentation  of  the  aerodynamic  parameters  as  well  as  saturation  limitations  tend 
to  support  the  conclusions  with  respect  to  the  superiority  of  the  proportional 
navigation  law.  It  was  found  in  this  study  that  the  use  of  a  "noise  adaptive" 
gain  in  conjunction  with  the  proportional  navigation  law  greatly  improved  the 
performance  of  the  systems.  The  noise  adaptive  gain  mechanization  was  ex¬ 
pressed  as 


_ G 

1  + 


v_ 

Vo 


.  .  o 

where  V  *  £  (A.  -  A) 

i-N-n 

o 

N  -  current  time  interval 


no  =  no  of  samples  in  window 

A  =*  ith  value  of  LOS  rate 
~  I  i-»  . 

A  =  —  £  A ,  average  value  of  the  LOS  rate  in  window 

i  no  .  „  i 
i=N-n 

o 

Vo  =  reference  variance 
G.  =  noise  adjusted  gain 

A 

G  51  nominal  gain 

For  periods  of  very  high  noise  levels,  this  law  was  quite  effective  in  reducing 
the  loss  in  energy  produced  by  noise  induced  maneuver  drag-  Since  the  present 
studies  show  that  a  reduction  in  gain  reduces  the  sensor  noise  contribution  to 
miss  distance,  this  should  also  provide  some  insulation  against  degradation 


caused  fay  sensor  noise.  Care  must  be  taken,  however,  that  a  proper  balance  is 

struck  since  a  reduction  in  gain  also  produces  an  increase  in  contribution  of 
gyro  drift. 

Since  the  noise  adaptive  gain  is  relatively  easy  to  implement  in  a  system, 
development  of  methods  for  optimizing  its  parameters  Vq  and  nQ  would  be  a  good 
direction  in  which  to  proceed. 
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